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Multiterminal measurements have typically been employed to examine electronic properties of 
strongly correlated electronic materials such as transition metal oxides without the influence of 
contact effects. In contrast, in this work we investigate the interface properties of Fe3C>4 with 
different metals, with the contact effects providing a window on the physics at work in the correlated 
oxide. Contact resistances are determined by means of four-terminal electrical measurements as a 
function of source voltage and temperature. Contact resistances vary systematically with the work 
function of the electrode metal, <j)(M), M =Cu, Au and Pt, with higher work function yielding lower 
contact resistance. This trend and the observation that contact resistances are directly proportional 
to the Fe304 resistivity are consistent with modeling the oxide as an effective p-type semiconductor 
with hopping transport. The jumps in contact resistance values at the bias-driven insulator-metal 
transition have a similar trend with <f>(M), consistent with the transition mechanism of charge gap 
closure by electric field. 

PACS numbers: 71.30.+h,73.50.-h,72.20.Ht 



Transition metal oxides (TMO) with strong 3d- 
electron correlations exhibit a rich variety of physical 
properties [l|, such as metal- insulator transitions, high 
temperature superconductivity, magnetic ordering, resis- 
tive switching, etc. Bulk electronic transport of such 
strongly correlated electronic systems (SCESs) have been 
examined extensively. However, there has been little 
study of the contacts between normal metals (M) and 
SCESs, even though injecting charge into a SCES in- 
volves "dressing" normal metal quasiparticles due to 
strong correlations. M/SCES interfaces under bias have 
been considered theoretically, and recent studies suggest 
SCESs with correlation-induced gaps in the density of 
states may be treated as classical semiconductors with 
minor modifications 0; correlation effects, however, are 
predicted to suppress rectification at M/SCES interface 
|3(. In such gapped SCESs, breakdown of the insulat- 
ing state upon application of an external electric field 
is theoretically described as analogous to Landau-Zcner 
tunneling in band insulators jj, which is effectively the 
closure of the charge gap by electric field Q. 

Contacts also play a critical role in resistive switch- 
ing (RS), in which a significant change of electrical re- 
sistance (from a comparatively insulating Off state to a 
more conductive On state) is induced by an application 
of external stimuli such as current, voltage or electric 
field @. Electric field driven breakdown of a gapped cor- 
related state has been proposed as the mechanism for 
RS in complex manganitespj and magnetite [§-[Toj . Non- 
volatile RS has also been observed in other TMOs, in- 
cluding perovskites[TIl, Ti0 2 [Hill, and doped SrTi0 3 
fl4j . Nonvolatile RS systems are of much technologi- 
cal interest fl5| . Experiments probing the M /TMO in- 



terface properties in nonvolatile RS systems are scarce. 
Multilead electrical measurements [16( and impedance 
spectroscopy [l?) studies show that the RS takes place 
at the interface. RS behavior is dependent on the work 
function of the metal in contact with the TMO [IH, [lj| ; 
RS is observed only if the contact metal forms a rectifying 
contact with the oxide material [(I |20|. 

We present a systematic experimental study of inter- 
face transport properties between a SCES TMO, mag- 
netite (E93O4), and three metals, M =Cu, Au, Pt, using 
nanostructure techniques. Stoichiometric bulk magnetite 
undergoes a phase transition between a high tempera- 
ture "bad metal" phase and a low temperature insulating 
phase at the Verwey temperature, Ty ~ 122 K[2l[, and 
exhibits electric field driven RS below Tv[3-[l0j. We use 
multiterminal measurements to assess contact resistances 
as a function of bias voltage, temperature, and work func- 
tion, </>, of the contacting metal [0(Cu)<0(Au)<</>(Pt)]. 
For all contacting metals, low bias contact resistances 
are found to track the bulk resistivity as a function of 
temperature, with larger <f> correlating with lower con- 
tact resistance. We also track contact resistances and 
their changes through the RS transition. These data are 
consistent with field-driven gap closure as the RS mech- 
anism. 

The 50 nm Fe304 (100) thin films used in the present 
study were grown on (100) oriented MgO single crystal 
substrates as described elsewhere [22|, [23| . Devices for 
four-terminal measurements were prepared on the sur- 
face of the Fes04 films by electron beam lithography and 
electron-beam deposition of contact metals, as described 
in detail in Ref. [l0( . To rule out sample-to-sample vari- 
ations, sets of devices were made incorporating differ- 



2 




FIG. 1. (a) NPGS pattern of different contact metal devices 
on the same piece of Fe3C>4 film. Grey squares are 300 /an 
x 300 (jm pads leading to source and drain electrodes, while 
pads colored in white lead to the voltage probes. The very 
small box in the center of the upper-left device defines a re- 
gion magnified in (b). Colored in black and grey, the [- and 
J -shaped source/ drain leads are the ones made of different 
contact metals, (c) and (d) are SEM images of four-terminal 
devices with Au-Au and Cu/Au-Pt source-drain leads, respec- 
tively. 

ent combinations of contact metals as source and drain 
electrodes: Cu-Au and Cu-Pt, as well as Au-Au and Pt- 
Pt, on a single magnetite wafer. Thus all devices under- 
went the same fabrication steps (Fig.QJ,), and all devices 
had the same geometrical characteristics. The Verwey 
temperature of the resulting Fe3 04 devices, found from 
zero-bias resistance, i?, as a function of temperature, T, 
is ^95 K, well below the Ty «110 K of pristine films 
grown under identical conditions. This difference is due 
to the intensive fabrication procedure, likely from devia- 
tions from the ideal FeaCu stoichiometry, known to sup- 
press Tviiimi. 

The whole four-terminal pattern except for the \- and 
J-shaped source/drain leads (see Fig. [lb), are made with 
a 6 nm Cu adhesion layer[l(| and a 10 nm Au cover layer. 
The [ and J source-drain leads are made of 20 nm layers 
of either Cu, Au or Pt (Fig. [lj>). This design ensures as 
much as possible that all devices have identical voltage 
probes, and that the difference in performance is indeed 
due to the difference in the interface with the source and 
drain contact metal. Some devices have shared large con- 
tact pads (Fig. [Th). 

Fig. QJi shows the device layout. Representative SEM 
images of Au-Au and Cu/Au-Pt devices are shown in 
Fig. [TJ; and d, respectively, demonstrating source and 
drain electrodes along with two pairs of voltage probes 
within the ^1/^m-long channel. Gold wires (25 fim di- 
ameter) were attached by In soldering to the 300 /im x 
300 /im pads of the devices to connect to the pins of the 
measurement puck. The puck was placed into a cryostat 
(300 K - 80 K) for measurement. Electrical characteriza- 
tion was performed via standard four-terminal methods 
using a semiconductor parameter analyzer (HP 4155 A). 
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FIG. 2. Schematics of (a) electrical circuit for four-probe 
measurements. Letters S and D denote source and drain con- 
tacts, respectively, and (b) how the source voltage, Vout, is 
distributed along the channel, (c) Temperature dependence 
of the low bias device resistance, -Rdev, and contact resistance 
Rc(M), M=Cu, Au, Pt, demonstrating R c (Pt)< Rc(Au)< 
iic(Cu), and that for all cases Rc ~ -Rdev- Inset shows the 
weak dependence of Rc (M) on source voltage in the (nearly) 
linear regime around V at 80 K. 

The voltage, V ou t, is applied to the source lead with the 
drain grounded, and current flowing through the channel 
is recorded, while a pair of voltage probes senses voltages 
V\ and V2 (and therefore AV = V\ — V2), as shown in 
Fig. Hi. 

The total applied potential difference, V^ut, must equal 
the sum of voltage drops across a device. There is 
a voltage across the source electrode/Fe3C>4 interface 
[Vc (source)] (Fig. [2b). In the assumption of a homo- 
geneous film exhibiting ohmic conduction between the 
electrodes, voltage linearly drops across the channel. Fi- 
nally, there is a further voltage drop at the FeaC^/drain 
electrode interface [Vc (drain)] ((Fig. HJd)). With these 
assumptions and using the geometrical characteristics of 
our devices as obtained from SEM images (i.e., £1, A£, 
and £2, see Fig. [Tb) we can calculate Vc (source) and 
Vc (drain) [l(| at each V out . 

Four measurements were conducted at each T for each 
four-terminal device with two pairs of voltage probes. 
With source and drain chosen, measurements were made 
with each pair of V probes; source and drain were then 
reversed, and the V probe four-terminal measurements 
were then repeated. From the resulting data, Rc (source) , 
-Rdev an d Rc (drain) (the source contact resistance, the 
device resistance, and the drain contact resistance, re- 
spectively) were calculated at each V ou t according to the 
procedure described below and in Ref. [Tcj. Values ob- 
tained with different source/drain and voltage probe con- 
figurations were checked for consistency. Calculations 
give the same values of i?DEV; -Rc (source) and Rc (drain) 
for left and right pairs of voltage probes, and Rc (source) 
and Rc (drain) switch appropriately upon exchanging the 
biasing/grounding choice. 

Wc first consider the properties of the M/Fe304 inter- 
face near zero bias. Around V (typically -1 V to 1 V 
range of V ou t) current-voltage (I-V ou t) characteristics of 
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the devices are linear. In this Ohmic regime the contact 
resistances at each V^ ut can be calculated as Rq (source) = 
Vc (source)// and i?c(drain)= Vc (drain)//, where / is a 
measured current at a certain V ou t- This data is shown 
as a function of temperature in Fig. [2J; for a particu- 
lar source/drain and voltage probe configuration, mea- 
sured with V^ u t | < 1 V. The Vcrwcy transition is ap- 
parent as the kink and change in slope at 95 K. The 
trend i?c(Pt)< -ffc(Au)< i?c(Cu) is consistent for all 
temperatures. It is striking that Rq(T) tracks -Rdev(T'), 
proportional to the magnetite resistivity, over the whole 
temperature range. We discuss this further below. 

The inset of Fig. |2J: demonstrates the weak voltage 
dependence at low bias of Rq for Cu, Au and Pt at 
80 K, extracted from the probing of Cu-Au and Cu- 
Pt devices. The plots clearly illustrate that Rc(Pt)< 
i?c(Au)< Rc(Gu). Note, that i? c (Cu) extracted from 
the independent measurements of Cu-Au (closed red 
dots) and Cu-Pt (open black squares) devices is essen- 
tially identical. 

The observed trend of Rq with metal type is consis- 
tent with the energetic offset between the Fermi level 
of the metal and the states responsible for transport in 
the magnetite. In the absence of interfacial dipoles and 
Fermi level pinning, this offset is determined by the dif- 
ference in the work function, fa of a contact metal M 
and FC3O4. The work function of clean magnetite is 
</Wfc(Fc304)=5.78 cV, while (f> 01 clean Cu, Au and Pt 
are 4.65 eV, 5.1 eV and Pt 5.65 cV, respectively (26} . 
While work functions can change due to surface ter- 
mination and adsorbates[27j, the systematic trend sug- 
gests that the ordering of work functions remains the 
same. Magnetite surfaces were freshly cleaned immedi- 
ately prior to insertion of the samples into the metal de- 
position system. For our metals, all of the M/FC3O4 
contacts appear to have 0(Fe3O4)> <fi(M). 

The dominant type of charge carriers in Fe304 de- 
pends on the stoichiometry of the material as well 
as the temperature range [28l [29j . In bulk stochimet- 
ric magnetite, below Ty the charge carriers (as deter- 
mined via Hall effect measurements) are holes. Fe3C>4 
might therefore be modeled as a p-type semiconductor 
(s/c)[28|,[3(|, though it is important to keep in mind that 
the charge carriers in magnetite are not well-described by 
the picture of extended states in conventional crystalline 
semiconductors |2 1| . The hopping character of charge 
transport in magnetite both above and below the Verwey 
transition temperature suggests that a more apt com- 
parison would be with highly disordered semiconductors, 
such as organic polymers (3li [32l ]. In such hopping sys- 
tems, when a carrier is injected from a metal electrode, 
there is a competition between diffusion of the carrier 
away from the contact by hopping, and the attractive 
interaction between the localized carrier and its image 
potential in the metal electrode. As a result, lower bulk 
resistivity leads to lower effective contact resistance [3lj. 

A schematic of the suggested interface energetics is 
shown in Fig. [3] The key difference between this and the 
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FIG. 3. Schematic diagram of the suggested energetics at the 
metal/oxide interface below the Verwey temperature. Various 
parameters are defined in the text. Charge transport takes 
place between the single-particle states at the metal Fermi 
level and localized states in the broadened "valence" band 
of the oxide. Above the Verwey temperature, the situation 
is similar, with transport in the oxide still taking place via 
the hopping of localized carriers. The trends in the contact 
resistance data suggest that in our samples dominant carriers 
are hole-like throughout the whole temperature range. 

conventional Schottky case is that the electronic states 
in the broadened "valence" and "conduction" bands of 
the oxide are localized, so that conduction takes place 
via hopping [2l|. In the framework of conventional bar- 
rier formation at M/p- type s/c interface with fap-tyj>c 
s/c)> 4>(M), the height of the Schottky barrier, fa is: 

fa = E g -[<f>(M)- X ]=E g + X -<j>(M) (1) 

where E g and \ are a band gap and an electron affinity 
of a semiconductor, respectively [33] . According to Eq. [TJ 
the higher (j>(M) the lower the barrier height, fa and, 
thus, the lower Rq of the contact interface. This is con- 
sistent with the observed trend in Fig.[5J: below Ty. The 
further observation that Rc{T) oc i?DEv(/ 1 ) is consistent 
with charge injection into other hopping systems (3ll. |32| . 
Note that this does not imply that there is no change 
in conduction across the Verwey transition. Rather, it 
implies that transport in the oxide both above and be- 
low the transition temperature involves hopping of charge 
carriers, rather than transport via extended states. 

It is important to note that the effective bulk resistivity 
of magnetite films may be strongly affected by antiphase 
boundaries (APBs), unavoidable at some density in films 
due to the kinetics of nucleation and growth. APBs are 
observed in Fe^O^ epitaxial films grown onto MgO (100) 
substrates [13, H^]. Since the Fe304 (Fd3m) crystal struc- 
ture is of lower symmetry than MgO (Fm3m) , there are 
several equivalent nucleation sites on the MgO (100) sur- 
face, which leads to the formation of APBs at the inter- 
face of neighboring grains. The presence of APBs leads 
to a large difference in the physical properties of epi- 
taxial magnetite films and those of the bulk, including 
an enhanced resistivity |35j. magnetoresistanceprjj . and 
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magnetization that does not saturate at high magnetic 
fields d?}- The complicated magnetic microstructures oc- 
curring due to the APB network is likely to reduce the 
spin polarization of magnetite films [38j. though spin po- 
larization is not relevant to the present work. In a care- 
ful study of the thickness dependence of magnetite film 
properties, Eerenstein et al. observed[35j that APBs con- 
tribute significantly to effective resistivity at thicknesses 
below 30 nm. In our 50 nm films, measured resistivity is 
comparable to bulk resistivity, suggesting that APBs are 
not the dominant source of scattering. It would be very 
interesting to perform contact resistance measurements 
in very thin films (in the limit that APBs strongly mod- 
ify the bulk resistivity), to see whether the trend shown 
in Fig. [2j; continues to hold; one would expect not. 
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FIG. 4. I-V ut (black) and I-AV (red) characteristics at 80 K 
(a) in non-linear regime before transition and (b) if applied 
voltage exceed a critical switching value, V sw , inducing resis- 
tive switching. 

Now we consider M/Fe3C>4 contacts at higher biases, 
far away from equilibrium. As the applied voltage in- 
creases, I-Vout curves become non-linear, but remain 
symmetrical around V (Fig. 0^) , in contrast to the rec- 
tifying curves of typical metal/semiconductor junctions. 
At the same time, I-AV curves remain linear (Fig. [4^,) 
even at high biases, suggesting that the non-linearity in 
I-Vout originates from the M/Fe304 interfaces. The ab- 
sence of significant rectification is consistent with expec- 
tations from strong electron-electron correlations and is 
typical of TMOs H. As V on t increases further, as re- 
ported previously [8j-ll0[ there is a sharp jump in current 
as V out reaches a critical switching value, V sw (Fig. 0b). 
The system after this RS transition is characterized by 
an approximately linear I-V ou t dependence and lower re- 
sistance compared to the state before transition. RS in 
magnetite devices is observed for contacts made from a 
variety of low- and high-work function metals (4.65 eV - 
5.65 eV range for Cu, Au and Pt, as well as Fe, Ti and Al, 
which have even lower 4>s of 4.5 eV, 4.33 eV and 4.28 eV 
[26j . respectively). 

The observation of reproducible field-driven RS in 
magnetite across a wide variety of contact metals is con- 
sistent with the proposed field-driven gap closure mech- 
anism, and stands in sharp contrast to trends seen in 
nonvolatile RS systems. In such systems, for which the 



origin of RS is attributed to the change in Schottky bar- 
rier height (or width) by trapped charge carriers at the 
interface states [HI, switching occurs only if the contact 
metal forms a rectifying contact with the oxide material. 
That is, p-type materials (such as complex manganitcs 
i?eo.7A .3Mn03 (Re = Pr, La; A = Ca, Sr)) demonstrate 
RS with low work-function metals, M=Mg, Al and Ti, 
and do not show RS with high work function metals such 
as Au and Pt f2Cj. On the other hand, n-type materi- 
als (Nb or Cr doped SrTiOs) demonstrate RS only when 
electrodes are made of high cf> material, such as Au or 
SrTiOa The RS in magnetite, in contrast, is clearly 
a bulk rather than a contact phenomenon. 

In the non-linear regime the definition of contact resis- 
tance as i?c = Vol I is no longer physically meaningful. 
In this case we estimate differential Rq as AVc / Al con- 
tinuously over narrow AV ou t ranges (~ 0.5 V) in which 
I-Vout curves are approximately linear. In the apparently 
linear regime at biases above V SV! , i?c is again calculated 
the same way as in linear regime around V. 
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FIG. 5. (a) Calculated differential Rc(M), M=Cu, Au 
and Pt, and resistance of the magnetite channel, -Rdev, 
at 82 K before and after transition point. (b) Temper- 
ature dependence of the absolute values of jumps in con- 
tact resistances, \R 3 Q mp (M)\. Error bars represent a stan- 
dard deviation over several independent measurements (dif- 
ferent devices, left/right voltage pairs, different grounds, pos- 
itive/negative switching branches). 

The representative dependences of calculated differen- 
tial Rq on V ou t for Cu, Au and Pt contact metals are 
shown in Fig. [5^ at 82 K (only positive voltage branch 
for clarity; negative branch is symmetrical). Reflect- 
ing the non-linearity in I-V ou t curves, contact resistances 
decrease as V out increases. The Rc(Pt)< i?c(Au)< 
i?c(Cu) dependence is observed in the whole V out range 
up to the transition point (Fig. [SJ. 

On the onset of the transition, as V ou t reaches criti- 
cal V sw value, contact resistances decrease abruptly and 
Rc values in the On state are approximately identical 
(about 600-700 SI for this set of devices) for all three 
contact metals (Fig. [5]). Note that contact resistances of 
both source and drain electrodes drop at the transition 
point. The resistance of the magnetite channel, i?DEVi 
also decreases at V out =V sw (Fig. [SJi). The absolute value 
of i?DEV> however, still scales with the channel length 



5 



(for details see Ref.[10(), as expected. It is clear from 
Fig. [S^, that i?DEV vs Kut dependence is "natter" (closer 
to Ohmic), than the contact resistance's. This is a direct 
manifestation of the linearity of I-AV curves in contrast 
to the non-linearity of I-V ou t curves (see Fig.|4^,). 

In previous work[t|[l(3| we proposed chargegap closure 
by electric field recently predicted in theory 0,[5[ as con- 
sistent with RS in magnetite. In this mechanism, with 
the assumption of Schottky-like barriers at M/Fe304 in- 
terfaces, the magnitudes of contact resistance jumps at 
a transition point, R° ( ^ mp , are expected to depend on 
the work function of the metal according to the relative 
alignment of Fermi level, Ep, of the contact metal and 
effective E F of Fc 3 4 . 

The data presented experimentally demonstrate 
that i?° // (Pt)< R% ff (Au)< R° If (Cu) at all V out 
before the transition, and i?g"(Pt)« i?g n (Au)« 
i?g"(Cu) after the transition (fig. Eh). Thus, 
jumps in R j g mp \M)=R% n {M)-R° } ' } \M) arc indeed in- 
dependent as Rg mp {Pt)< R^ mp {Au)< R% mp (Cu). 
However, as can be seen in Fig. [5^, the value of R 3 ( ^ mp 
is also influenced by the V sw value which, in turn, is pro- 
portional to the channel length and thus varies from de- 
vice to device. That is the reason that R 3 Q mp values 
vary significantly (see error bars in Fig. [5Jd) in different 
measurements of different devices. In some cases this 
variation obscures the trend R J ( ^ mp {Pt)< R J Q mp (Au)< 
J?p"" p (Cu), but overall the latter dependence stands true 



and consistent with expectations of the charge gap clo- 
sure mechanism. 

We have investigated M/Fe 3 4 (M=Cu, Au and Pt) 
contacts in this correlated resistive switching system by 
means of four-terminal electrical measurements. The 
dependence of contact resistances on metal work func- 
tion <j> is consistent with modeling the correlated oxide 
as an effective p-type semiconductor. Likewise, the ob- 
served proportionality between contact resistances and 
bulk resistivity over a broad temperature range is con- 
sistent with models and observations of charge injection 
into semiconducting hopping systems. We continuously 
traced the contact resistances as a function of bias before 
and after the resistive switching transition and showed 
that the jumps in contact resistances evolve with (f> as 
4™ p (Pt)< Rg mp {Au)< Rg mp (Cu) at all tempera- 
tures below Ty when switching is observed. This trend is 
consistent with the charge gap closure by electric field as 
a proposed driving mechanism of RS in magnetite. These 
results demonstrate the utility of nanostructure contact 
studies in unraveling the physics at work in correlated 
materials. 
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